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ABSTRACT 



Aims. By detecting high frequency quasi-periodic oscillations (QPOs) and estimating frequencies of them during the decaying tail of 
giant flares from Soft Gamma-ray Repeaters (SGRs) useful constraints for the equation of state (EoS) of superdense matter may be 
obtained via comparison with theoretical predictions of eigenfrequencies. 

Methods. We used the data collected by the Rossi X-Ray Timing Explorer (RXTE/XTE) Proportional Counter Array (PCA) of a giant 
flare of SGR 1806-20 on 2004 Dec 27 and applied a Bayesian periodicity detection method (Gregory & Loredo, 1992) for the search 
of oscillations of transient nature. 

Results. In addition to the already detected frequencies, we found a few new frequencies (f^pos ~ 16.9, 21.4, 36.4, 59.0, 1 16.3Hz) of 
oscillations predicted by Colaiuda et al. (2009) based on the APR 14 EoS (Akmal et al., 1998) for SGR 1806-20. 
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1. Introduction 

The study of periods of activity of Soft Gamma-r ay Repeaters 
(SGRs, for a recent review see [Mereghetti 2008) showing re- 
current bursts with sub-second duration and much more extreme 
events known as giant flares (on rare occasions), may have an 
important input to our understanding of neutron stars (NSs). 

In particular, the detection and analysis of the pres- 
ence of quasiperiodic oscillations (QPOs), up to a few 
kHz, has triggered a number of theoretical studies for pre- 
diction and direct comparison with various equations of 
state for superdense matter. These oscillations have been 
interpreted initially as torsional oscillations o f the crust 
dSamuelsson & Anderssonl2007l;ISotani et al.ll2007l) . While later 
the Alfv en oscillations of the fluid core have been taken into 
account (lLevinll2007l; ISotani et al.l[2008bt IColaiuda et al.ll 2009 ; 
Cerda-Duran et al. 2009; van Hov en & Levinl2010tlGabler et al.l 
2010; IColaiuda & Kokko tas 2010). These studies suggested how 
the observations can constrain both the mass, the radius, the 
thickness of the crust and the strength of magnetic field of NSs. 
In particular, the timing analysis of the decaying tail of the un- 
precedent giant flare of SGR 1806-20 on 2004 Dec 27, allowed 
to detect QPO frequ encies approximately at 18, 26, 30, 92, 150 , 
625, and 1840 Hz (llsrael et al Jl2005HWatts & Strohmaverl 2006: 
IStrohmaver & Wattsll2006HWatts & Strohmaverll2007l) in differ- 
ent time intervals, different rotation phase and different ampli- 
tudes of oscillations, by means of computation and analysis of 
the averaged power spectrum. 

Here, we present the results of a Bayesian approach of timing 
analysis, of the giant flare data set of SGR 1806-20 registered 
on 2004 Dec 27 by RXTE PCA, for detecti ng a periodic signal 
of unk nown shape and period developed by iGregorv & Loredol 
( 19921) (for its sensitivit y and advantages see further and also 
Grego ry & Loredoll 1996ft . 



2. Observational data and results of analysis 

The giant flare of SGR 1806-20 has been observed by many 
space based missions. The data recorded by RXTE PCA instru- 
ment, consisting of five Xenon-filled detectors covering the en- 
ergy range 2-50 keV, used the configuration GoodXenon, which 
records all good events detected in the Xenon chamber with full 
timing accuracy of 1 [is. Publicly available data were retrieved 
by th e XTE Data Finder (XDF) user interface dRots & Hilldrupl 
1997). Event data files have been created and photon arrival 
times were corrected for the solar barycenter using scripts pro- 
vided in the package XTEftools. The data set consists of 698770 
registered photons. However, during the initial intensive spike 
phase of the giant flare the detectors were saturated. For that 
reason, in our analysis we use the data after ~ 8.9sec of the flare 
onset, consisting of ~ 650000 photon arrival times, clearly cov- 
ering 51 rotational cycles of SGR 1806-20 (see Fig. |TJ. 

It is clear that observational detection and parameter estima- 
tion of QPO frequencies may play a crucial role for testing any 
theory predicting eigenfrequencies of the neutron star. In this 
connection, timing analysis of complex flare data set of SGR 
1806-20, with the aim of QPO detection, may be divided into 
several mutually connected, challenging problems. They include 
the significance of quasi-periodic signal detection and parameter 
estimation with high precision. Indeed, the decaying tail, of the 
giant flare of SGR 1806-20, itself has a bumpy structure, very 
complex light curve shape modulated by rotation of the neutron 
star (see Fig.fTJ. 

The most widely-used procedure for detection of QPOs is an 
analysis of the power spectrum calculated from the Fast Fourier 
transform (FFT) of uniformly sampled data. Application of the 
FFT to unevenly spaced arrival time series data requires binning 
of the data to produce equally spaced samples. Binning is a sub- 
jective procedure; the choice of the bin width and edges can af- 
fect the apparent significance of a detection and limits sensitiv- 
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tail of the flare (iBretfhorstll 1 9881: iGregorv & Loredol 1 992l [l 996: 



Fig. 1. Light curve of the RXTE giant flare of SGR 1806-20 ob- 
served on 27th of December 2004. The general decay of the gi- 
ant flare with a bumpy structure on top of it, namely a strongly 
periodic signal due to rotation (7.56 sec). The inset panel shows 
the rotational modulated light curve (filled circles) toget her with 
fitted piecewise constant model (solid line, Hutter 2007), shown 
here only for 2.5 rotational cycles, with a very complex light 
curve structure (see text for details). 



ity on short time scales. Moreover, the presence of the so-called 
red-noise (low frequency signal in the data, i.e. large calculated 
power can appear at low frequency because of long-timescale 
features of the data and at high frequencies from higher harmon- 
ics of complex shape strong periodic signal), may ca use certain 
probl ems with the interpretation of the results (e.g. iBretthorsil 
[l988h . In the previo us timing analysis of t he giant flare data set 
of SGR 1 806-20 bv llsrael et al] d2005l) and lWatts & Strohmaverl 
(120061) . an averaged power spectrum was considered. Namely, 
for each rotational cycle or certain rotational phase interval of 
SGR 1806-20 an independent classical power spectrum was de- 
termined depending on the phase of rotation and decaying tail 
of the giant flare, which were co-added and averaged subse- 
quently. This approach divides the data set into small time in- 
tervals, which automatically reduces the significance of any pe- 
riodicity detection of transient nature and may cause problems 
related to the wind owing of the data sets (Ivan der Klisl ll989: 
iGalleani et alJ [2001). The complex shape of the light curve (ro- 
tational modulation and decaying tail of flare light curves, Figjl) 
is not taken into account in a proper way, i.e. assuming that dur- 
ing small time/phase intervals there are no significant intrinsic 
variations in the data subsets. Moreover, in some circ umstances 
FFT may fail to detect the p eriodic signal (see e.g. Bret thorsll 
988; Gregory & Loredoll 19961). As explicitly shown by I Javnesl 
( 19871) (see also. lBretthorstl'l988ll2001l:lGregorvl2005h the prob- 
ability for the frequency of a periodic sinusoidal signal is given 
approximately by 



p(f n \DJ) oc exp 



C(f„) 



cr- 



(1) 



where C(f n ) is the squared magnitude of the FFT, show- 
ing that proper approach to convert C(/„) into probabilities in- 
volves first dividing by the noise variance and then exponentiat- 
ing, which suppresses spurious ripples in the power spectrum. 

We used a procedure which does not require binning and 
takes into consideration the rotational modulation and decaying 



iJavnes & B retthors i2003UVaughanl2009l) . 

For the analysis of the data for the search of QPOs during 
the giant flare of SGR 1806-20, we app lied a Bayesian method 
developed by Greg ory & Loredol (H 992) (hereafter referred to as 
the GL method) for the search of pulsed emission from pulsars 
in X-ray data, consisting of the arrival times of events, when we 
have no specific prior information about the shape of the signal. 
As a particular case of QPO can be considered a periodic signal 
with some length of coherence (Q = vo/2cr, vo is the centroid of 
the frequency and cr half width at half maximum), i.e. a periodic 
signal with additional parameters of the oscillation with start and 
end times 0. 

The GL method for timing analysis first tests if a constant, 
variable or periodic signal is present in a data set. 

In the GL method, periodic models are represented by a 
signal folded into trial frequency with a light curve shape as 
a stepwise function with m phase bins per period plus a noise 
contribution. With such a model we are able to approximate 
a phase folded light curve of any shape. Hypotheses for de- 
tecting periodic signals represents a class of stepwise periodic 
models with the following parameters: trial period, phase, noise 
parameter and number of bins (m). The most probable model 
parameters are then estimated by marginalization of the poste- 
rior probability over the prior specified range of each parame- 
ter. In Bayesian statistics posterior probability contains a term 
that penalizes complex models (unless there is no significant ev- 
idence to support that hypothesis), hence we calculate the pos- 
terior probability by marginalizing over a range of models, cor- 
responding to a prior range of number of phase bins, m, from 
2 to 12. Moreover, the GL method is well suited for variability 
detection, i.e. to characterize an a rbitrary sh ape light curve with 
piecewise constant function Z(f) (Rots 2006). 

For the search and detection of QPOs we used a slightly dif- 
ferent version of the GL method: First we determine Z(f) - fitting 
with a pie cewise consta nt model (solid line in the inset panel of 
theFig.m lHutterll2007h 

to characterize the complex light curve shape in the data 
set, giant flare decaying tail and rotational modulated light 
curve. Then we subsequently compare competing hypothesis, 
i.e. whether the data support a purely piecewise constant or 
piecewise constant +periodic model. If there is an indication of 
the presence of a periodic si gnal (e.g. odds ratio of competi ng 
models exceeding 1, see also [Gregory & Loredol (1 19921 1 19961) ^ 

we determine also the time intervals where it has its maxi- 
mum strength (e.g. amplitude or pulsed fraction) via a Markov 
Chain Monte Carlo (MCMC) approach using QPO start and end 
times as free parameters. 

Finally, in the latter case, we estimate parameters (frequency, 
phase, amplitude, coherence length of QPO, etc.) of the periodic 
signal with high precision. For example, in order to estimate the 
frequency of a periodic signal the posterior probabilities density 
function used: 
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p(a)\D,Mm) = - . ... . 
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where C = 



Ja>i„ u JO 



and W m (co, <p) = 



^0 W m (a>,0) 

is the number of ways the binned distribution could have arisen 
"by chance", n ; (w, 4>) is the number of events falling into the j th 
of m phase bins given the frequency, lo, and phase, <f>, and N is 



More strictly: coherent transient signal dvan der Klisll2010h 
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Fig. 2. Application of the GL method for blind periodicity search 
to the complete data set obtained by RXTE PCA of giant flare of 
SGR 1806-20 2004 Dec 27 revealed a strong coherent signal 
at the frequency of 0.13219244 Hz and higher harmonics up 
to 100Hz. Bayesian posterior probability density vs frequencies 
in the range of 12.0-160 Hz is shown. The inset panel shows a 
zoomed part of it around 16.88 Hz. Arrows are indicating higher 
harmonics (from 125 to 132) of the fundamental frequency and 
the dash-dotted vertical line shows one of the detected QPO fre- 
quency in a short time interval (see text, for details). 



Time interval 180 — 400 sec (flare onset time = 0) 




60 80 100 
Frequncy [Hz] 

Fig. 3. The Bayesian posterior probability density vs. trial fre- 
quency of the periodic signal for the giant flare data set of 
SGR 1806-20. QPO frequencies already detected by averaged 
power spectrum a n alysis and also with other mission RHESSI 
dlsrael et all 120051; IWatts & Strohmaverl 120061) are marked as 
small lines at the top axis. Those QPO frequencies are also de- 
tected by us. In addition, we have detected several more QPO 
frequencies (21, 59 and 11 6 Hz) with the Bayes ian method, 
which were also predicted bv lColaiuda et al.l d2009l) . 



the to tal number of photons (for details, see iGregorv & Loredol 
[1991 . 

First, we applied th e GL method as implemented by 
iGregorv & Loredol (1 19921) to the 2004 Dec 27 giant flare whole 
data set of SGR 1806-20 observed by RXTE PCA and started 
the timing analysis by performing a blind periodicity search in 
the range of 12.0-160 Hz. Naturally, we found a very strong co- 
herent signal at ~ 7.56s, the pulsation period of the NS, followed 
by higher harmonics up to the 100Hz (see Fig. [2J. 

Next, we divided the data set into 5 1 rotational cycles (see 
FigJTJ. Each of these subsamples were treated as independent 



data sets. We determined the Bayesian probability densities ver- 
sus trial frequency. Final probabilities were derived in two ways; 
First, as the multiplication (the likelihood) of those independent 
Bayesian posterior probability density functions; and the sec- 
ond one simply summing those independent probabilities (see 
Fig®. 

This approach also revealed a number of rotational cycles 
within which the probability of the model of periodic signal is 
significantly higher than the constant one. Namely, during time 
intervals of 183.8-191.2, 244.3-251.9 and 259.4-267.1 sec- 
onds, from the flare onset, odds ratios of periodic vs piece-wise 
constant models are ~ 10, ~ 30 and ~ 197, correspondingly. 

In order to detect start and end times oscillations we in- 
cluded also observational data of neighboring rotational cycles 
(where an oscillation with that frequency was not detected) and 
performed the periodicity search for an expanded time interval 
with additional two free parameters with the MCMC approach 
with Metropolis-Hastings algorithm. As initial values of these 
(tQpostan and tQPOend) parameters served start and end times of 
an observation, satisfying the condition: tobs. start ^ tQPOstart < 
tQPOend ^ tobs.end- Starting with the abovementioned initial val- 
ues can be considered as a good strategy, since the detected sig- 
nal has a higher significance in a subinterval of the considered 
time interval and the fast convergence of the MCMC procedure 
already provided. This analysis via MCMC revealed even shorter 
time intervals within which periodic signal is stronger. The es- 
timates of oscillations frequencies and the corresponding 68% 
interval of the highest probability densities are presented in the 
Table[T](see, also Fig. EE). 

As shown by our intensive simulations, this approach readily 
detects time intervals where periodic signal has highest and low- 
est strength. Nevertheless, we performed also time-frequency 
analysis (see Fig. [5] and [6j as well as dynamic power spectrum, 
Ivan der Klisl IT989) and candidate QPOs fitted to a Lorentzian 
function. We estimated frequency centroids and full widths at 
half maximums. For significant detections they are presented in 
TableQ] The coherence values Q = vq/2<t ~60-200 are indicat- 
ing on the transient nature of almost perfect periodic signals. 

For those significant detections we have also fitted folded 
light curves with a template function, t[r(t) = YIt=\ Ajsin[iu)o(t - 
i/fi)], truncating the series at the highest harmonic that was statis- 
tically significant after performing an F-test (see lDaU'Osso et aT] 
I2003L for details). It turned out that oscillations at 16.9Hz and 
36.8 Hz adding the second harmonic does not improve the fit 
of phase folded light curve (signifcance levels, accordinglly are 
0.71 and 0.16), while in the case of 21.3Hz it significantly im- 
proves it (significance level of 0.0006. Thus, oscillations at the 
frequencies of 16.9Hz and 36.8Hz can be considered of sinu- 
soidal shape, and at 21.3Hz not. We obtained similar results by 
analysing also power spectra, i.e. fitting spectral profiles of them 
with sine function. 



2 Addition rule of probabilities: probability that the QPOs at the given 
frequency are present A?i or At2 or both time intervals, while likelihood 
(i.e. summed power spectrum) defines the probability of presence QPOs 
during Afi and At2 (i.e. multiplication rule of probabilities, see also, 
Eq.[3 

3 A time interval where odds ratio of the models of a periodic 
and piecewise constant models just exceeds a theoretical minimum, 
Odds ratio = 1, i.e. both models have equal probability to be true. 
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Detected QPO freq u encies not reported i n the 
12005b IWatts & S trohmaver 2006; 



Table 1. 

literature (llsrael et al.. 
IStrohmaver & Wattsll200 



fQPO [Hz] 


Time intervals 


Lorentzian 


FWHM[H 


( 68% credible region) 


of QPOs" 


centroids(Hz) 




16.88" (16.87- 16.90) 


259.4-267.1 


16.90±0.004 


0.12 


21.36" (21.35 -21.38) 


244.3-251.9 


21.34+0.003 


0.35 


36.84" (36.83 - 36.88) 


183.8-191.2 


36.88±0.004 


0.24 


59.04 (58.58 -59.28) 


146.0-176.2 






61.26 (61.25 -61.27) 


251.9-395.6 






116.27 (116.24- 116.28) 


168.7-198.9 







a Giant Flare onset time is set to 

b Highly significant detection (for details, see text) 




16.80 16.85 16.90 16.95 17.00 0.0 0.5 1.0 1.5 2.0 

Frequency [Hz] Phase 

Fig. 4. Bayesian posterior probability density vs frequencies in 
the time interval of 259-267 sec from the giant flare onset. By 
dashed vertical lines are indicated the 68% region of the high- 
est probability density. Odds-ratio of periodic vs constant model 
is ~ 200. The right panel depicts the phase folded light curve 
with frequency /qpo = 16.88Hz,it has almost perfect sinusoidal 
shape (for details, see text). 



3. Discussion 

We report the detection of oscillations of transient nature from 
SGR 1806-20 giant flare decaying tail recorded by the RXTE 
PCA by applying the GL method for periodicity search. We have 
confirmed the detections of QPOs at frequencies (in t he range of 
12 - 160Hz) reported by I Watts & Strohmaverl d2006l) and, in ad- 
dition, we found some more QPOs at fopo s = 16.9, 21 .4, 36.8Hz 
with corresponding odds ratios ~ 197, ~ 30 and ~ 10, in shorter 
time scales, i.e. within individual rotational cycles (see, Fig|7]i. 

These odds ratios, describing the significance of presence of 
a periodic signal, are sensitive to the frequency search range. In 
contrast, detected frequencies of QPOs, found by locating max- 
imums in the posterior probability density function, are insensi- 
tive to the prior search range of frequencies. We computed the 
uncertainty of /qpos at 68% confidence level by using this pos- 
terior probability density function. In addition to that, we have 
also estimated the significance of our QPO detection by an em- 
pirically determined cumulative distribution (see, FigJTJ. 

The great variety of neutron star oscillation modes 
dMcDermott et al.l 11985) is associated with their global struc- 
ture as well as their internal characteristics. Neutron star seis- 
mology has already been proposed as a tool to understand their 
internal structure ( Kokkotas et al ] 12001b . The various types of 



Fig. 5. Bayesian posterior probability density vs frequencies and 
the time of the giant flare onset. Clearly seen significant detec- 
tions of oscillations at 16.9,21.4 and 36.4Hz. Large Q = vo/2<x 
values (see, Table [TJ indicating on the coherent nature of tran- 
sient oscillations (for details, see text). 



Dynamic power spectrum 



DFT 

I 



Fig. 6. Dynamic power spectrum of giant flare of SGR1806-20 
in the frequency range of 16.5-17.5 Hz and in the time interval 
of 250-275sec (see, TableQ}. 



QPOs in cycles: 24 32 34 




Log(Odds ratio) 

Fig. 7. Empirical cumulative distribution function as a function 
of Odds ratio, i.e. the probability of periodic model vs constant 
one, on the base of simulations of QPO frequencies with differ- 
ent amplitudes and noise, as much as possible to the observed 
values in terms of number of registered photons, spanned times, 
etc.. Detected QPO frequencies indicated by vertical dashed 
lines, showing a high significance of detections. 



4 



V. Hambaryan et al.: Bayesian timing of giant flare of SGR 1806-20 



oscillations carry information about the equation of state, the 
thick ness of the crust, the mas s, radius and even the rotation 
rate (Gaertig & Kokkotasll201fJh . Still, modeling a truly realis- 
tic oscillating neutron star is difficult, however the potential re- 
ward is considerable. This has been demonstrated from the re- 
cent theoretical results for the QPOs which have been inter- 
preted as magnetoelastic oscillations. These calculations have 
shown how the observations constrain both the mass, the radius, 
the thickne ss of the crust and the strength of magnetic field of 
these stars dSqtani et al.ll2Q07t ISamuelsson & Andersson1l2007t 
IColaiuda et alJl2009HCerda-Duran et al.ll2009l) . While one could 
even set sever e constrains in the g eometry of the interior mag- 
netic field (see lSotani et alJl2008al) . 

In particular, by studying Alfven oscillations, in the stellar 
interior, in the absence of crust, Colai uda et al.l £2009) cou ld re- 
veal t he existence of Alfven continua (as suggested by iLevinl 
120071) . The edges of this continuum have been used to explain 
the observed QPOs. It is worth noticing, that QPO frequencie s 
found in this work have bee n predicted bv | Colaiuda et al.l J2009). 
More recent calculations by Colaiuda & Kokkotai ( 120101) taking 
into account the presence of crust, i.e. the coupled crust-core os- 
cillations, verified the presence of the continuum but in addition 
they revealed the presence of discrete Alfven oscillations as well 
as some extra discrete modes which can be associated with the 
oscillations of the crust. These new results explain the previous 
published QPO frequencies but also the ones presented here and 
set quite severe constraints in the parameters of the neutron star. 

4. Conclusions 

1. We found new oscillationsof transient nature applying 
Bayesian timing analysis method of the decaying tail of 
the giant flare of SGR 1806-20 2004 Dec 27, observed by 
RXTE PCA, not yet reported in the literature. 

2. Some of those oscillations frequencies (fQpos ~ 
17, 22, 37, 56, 1 12Hz) are predicted by the theoretical 
study of torsio nal Alfven osc illations of magnetars (see, 
Table 2, by IColaiuda et al.1 l2009h . suggesting APR X4 
(lAkmal et al.il 19981) EoS: 4 ! of SGR 1806-20. 

3. These preliminary results are very promising and we plan 
to extend our high frequency oscillations research (both the 
theoretical predictions as well as the observations) to both 
activity periods, as well as to the quiescent state of SGRs 
and AXPs, as well as to the isolated neutron stars with com- 
paratively smaller magnetic fields. 
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